Dynamics after Interaction Quenches in One-Dimensional Fermionic Systems 
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We show that the dynamics of quenches in one dimension far off equilibrium can be described 
by power laws, but with exponents differing from the fully renormalized ones at lowest energies. 
Instead they depend on the initial state and its excitation energy. Furthermore, we found that 
for quenches to strong interactions unexpected similarities between systems in one and in infinite 
dimensions occur, indicating the dominance of local processes. 
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In the last years, seminal experimental setups have 
been developed combining a very good decoupling of 
the quantum systems from their environment with a 
high degree of controllability of the system's parameters. 
This renders the observation of the temporal evolution of 
closed quantum systems for long times possible. Chang- 
ing the internal parameters provides tools to investigate 
systems out of equilibrium. In optical lattices a change 
in the trapping potential and the external frequency is 
used to manipulate the internal parameters such as hop- 
ping and particle- particle interaction [ll-Q . Furthermore, 
various pumb-probe experiments based on ultrafast spec- 
troscopy [H[ have been developed 

These developments have triggered extensive theoret- 
ical studies of physics far from equilibrium, based on 
a large variety of analytical and numerical tools 0, 0- 
Ull . Hil - tlit . The goal is to qualitatively understand and 
to quantitatively describe the dynamics of quantum sys- 
tems far from thermal equilibrium. How do such states 
evolve? Do they relax towards equilibrium? How does 
this happen and on which time scales? 

These issues are relevant in all dimensions. But so far 
the one-dimensional (ID) case has attracted the greatest 
interest. On the one hand, this is because the ID case is 
amenable to a larger variety of approaches. On the other 
hand, there are fascinating conceptual issues inherent to 
ID systems: (i) All extended non-trivial integrable quan- 
tum systems are ID and the existence of an exhaustive 
number of conserved quantities influences the dynamics 
strongly [3]. (ii) Gapless ID systems are successfully 
described in the low-energy sector by free bosonic field 
theories 

[IB-Ill- 

Thus the question arises whether the 
same or similar bosonic field theories are also able to de- 
scribe the nonequilibrium dynamics [23H24I ]. This point is 
related to (i) because a free bosonic theory is integrable. 
Still, relaxation in local quantities due to dephasing may 
occur 
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In this work we provide evidence that gapless bosonic 
field theories can indeed describe nonequilibrium dynam- 
ics in infinite ID systems implying power laws. For the 
accessible times, the appropriate exponents are not the 
equilibrium exponents if the system is strongly quenched. 



The effect goes beyond the substitution 7 — > 2j(l + 7) 
which appears in quenches of Luttinger li quid s (LL) 
where 7 is the equilibrium exponent [3, H3L l26l|. Fur- 
ther, we argue that the power laws observed at short and 
moderate times remain meaningful for larger times if the 
macroscopic excitation energy AE/L (L: system size) in- 
duced by the quench is sufficiently large because AE/L 
cuts off the renormalization (RG) flow which would yield 
the equilibrium exponents at its fixed point. This is our 
first key result. Our second key result ensues from the 
comparison of the ID dynamics to the one at infinite di- 
mension where dynamic mean field theory (DMFT) is 
exact 0, H, [ll| . Much to our surprise strong qualitative 
similarities occur. 

Quenches are an efficient way to realize states far from 
equilibrium [25j . Often, the interaction value is changed 
abruptly increasing it @, 0-11, 13, 27 or decreasing 
it [15(. We focus on interactions which are suddenly 
switched on and consider Fermi seas as initial states. 

Two generic fermionic lattice models are studied. The 
first are spinless fermions 
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with nearest- neighbor repulsion (NN), where cj (dj) cre- 
ate (annihilate) a particle at site j and hj = c}-c-. The 
second is the Hubbard model comprising spin a 4. 



- J E 



fct c- 
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and local repulsion. We study U(t) = Q(t)U > and de- 
fine the band width W — 4 J as energy scale. Below some 
quantities with spin a are denoted; for the corresponding 
ones without spin a is to be omitted. 

The approach used is an expansion of the Heisenberg 
equations of motion for cj a in powers of the time t [l3[ . 

By commuting H after the quench recursively with a 
we obtain differential equations for the prcfactors of the 
expansion of ct a (t) in monomials of c| CT and at t = 0. 
The application of one commutation is a 'loop'; up to 
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11 loops are realized. A conceptual asset is that the ap- 
proach works for the infinite lattice by exploiting trans- 
lational invariance, see Supplemental Material. 
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FIG. 1: (color online) Upper panel: Jump An(t) for the half- 
filled Hubbard model for various loop numbers. Lower panel: 
An(t) for increasing U (from top to bottom at small t) in 11 
loops. 

We focus on the momentum distribution rik[t) '■= 
(c\. a c k a )(t) where k is the wave vector fl3| . In particu- 
lar we study the jump An(t) :— nk F+0 (t) — nf CF _ {t) at 
the Fermi vector k-p. In the initial Fermi sea An(0) = 1 
holds. After the quench, An shows slow relaxation to 
zero @,[l3|,[27| and coherent oscillations [HUff. The up- 
per panel of Fig. [T] depicts the jumps for increasing num- 
ber of loops for the half-filled Hubbard model. Good con- 
vergence is obtained for 11 loops up to about t ps 10/ VF. 
The precise value up to which the data is reliable de- 
pends on the details. The lower panel of Fig. Q] displays 
the data for various values of U . 

In Fig. [2] the jump An(t) is shown for the spinless 
NN model. The data agrees very well with the data by 
Karrasch et al. in Fig. 1 of Ref. |23| obtained by time- 
dependent infinite-size DMRG. We compare our data to 
the analytic power law (l3j 



An(t) = 



{2v p tf 



7p(1+7p) 



{2 V<7 tf 



7<t(1+7„) 



(3) 

where v 6 {p, a} stands for charge (p) or spin (er) chan- 
nel. The parameters r„ are length scales, v v velocities 
and 7^ the exponents which read 7„ = [K v + K" 1 — 2)/4 
as function of the K,, value from bosonization. For the 



spinless case all quantities coincide K = K p 



K„ 



(7 



and so on. The dashed lines in Fig. [2] use 



K GS = 7r/[2(7r - arccos(2E//W0)] (4a) 
v GS = 7rsin(arccos(2J7/W))/[2arccos(2{7/W)], (4b) 
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FIG. 2: (color online) Solid lines: An(t) in the half-filled 
NN model for increasing U (from top to bottom) in 11 loops. 
Dashed lines: An(t) as in Eq. (|3)) with the ground state (GS) 
exponents from Eq. @. Dashed-dotted lines: An(t) with the 
Fermi sea (FS) exponents from Eq. (JSJ). 



from Bethe ansatz [17|,[28[ reflecting the equilibrium exci- 
tations close to the ground state (GS). The cutoff length 
r is fitted and evolves from 0.4 and 0.6 on increasing U . 
The formulae are reasonable only up to U = W/2 where 
the system enters a gapped phase. 

The dashed-dotted lines use the values [22[ 



Kfs = y/(nvF-U)/{-KVF + 3U) 

V FS = (l/w)y/(TTV F +U) 2 -4U 2 



(5a) 
(5b) 



which result from bosonization around the Fermi sea (FS) 
IE El US [13 . The fitted cutoff length r varies from 0.4 
and 0.2 on rising U . The FS parameters characterize the 
LL Hamiltonian which describes bosonic modes in the 
vicinity of the FS, if processes not captured by the LL, 
e.g., the curvature of the dispersion, do not matter much. 
Then the FS exponents are surely relevant at short and 
moderate times when the system is still close to the FS. 
The GS parameters characterize the LL Hamiltonian of 
the bosonic modes close to the GS. Hence they are surely 
relevant for long times in equilibrium. 

The FS formulae are reasonable only up to U « W 
where «fs vanishes at U = irvp- For large U the effects 
of umklapp scattering and dispersion curvature hamper 
the bosonization approach. 

Inspecting Fig. [5] we see that the microscopic model 
displays oscillations in An(t) which are absent in Eq. 

These oscillations can be ascribed to the momentum 
cutoff of the interaction in microscopic models 0, [2(| . 
Otherwise, the power law © nicely describes the dy- 
namics at moderate times in agreement with Ref. I23I . It 
cannot be decided which of the two sets (GS and FS) 
fits better to the microscopic data, though the FS power 
law fits slightly better at U > 0AW. The reasons are 
the small differences in K from @ and ([S]) as long as 
U < W/2. This is accord with the closeness of the 
quenched NN system to equilibrium as measured by en- 
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FIG. 3: (color online) Left panel: Energy expectation value of 
the Fermi sea relative to the ground state energy at half-filling 
for the Hubbard and the NN model. Right panel: Period of 
coherent oscillations in the strong interaction regime of the 
Hubbard model; n is the filling factor of one spin species. 



ergy AE = (FS\H(t > 0)|FS) - E a which is excited in 
the system by the quench, see left panel of Fig. 31 The 
ground state energy E results from Bethe ansatz 28l.l29l]. 
For the NN model AE stays very small up to U = W/2 
where the gapped phase begins [28j| . For the Hubbard 
model much larger energies are reached. 
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FIG. 4: (color online) Jump An(i) for the Hubbard model at 
quarter-filling at U = 0.8W. Lines as in Fig. H For the GS 
exponents r p ~ 0.01 is fitted; r a does not occur. For the FS 
exponents r p ~ 1 and r a w 0.6 are fitted; 



To distinguish between GS and FS exponents on the 
available time scales a model is needed where they differ 
significantly. The Hubbard model is a good candidate 
because of its AE values (Fig. |3j and because one of its 
if -values (K a for the spin channel) differs sizeably for 
GS (K a ,GS = 1) and FS [H[2l| 



Kp,FS 

V P ,FS 
W<r,FS 



v / 2ttv f /(2ttv f + 2U) (6a) 

^2ttv f /(2ttv f - 2U) (6b) 

v FV /l + U/(nv F ) (6c) 

vpy/1 - U/(nv F ). (6d) 



On the time scales studied, the FS power law captures 
the microscopic dynamics of An(t) much better than the 
GS power law, see Supplemental Material for other U 
values. This supports our argument that the dynamics 
at short and moderate times can be described by a LL 
model with exponents relevant for the vicinity of the ini- 
tial state, here the FS. This is the first key result. 



Only the oscillations due to momentum cutoff j23l |26| 
are missed by the power law ([3j). For small values 
U Q.3W the differences between the GS and the FS 
power laws become small so that they are indistinguish- 
able. For larger U 2; W the agreement between mi- 
croscopic data and FS power laws deteriorates due the 
neglected physical processes (dispersion curvature, back 
scattering) in the LL. 

The next intriguing question is what happens for 
longer times tl To date this question cannot be answered 
definitely. Analytically, no RG theory for nonequilibrium 
physics exists. Numerically, no tools can treat long times 
far off equilibrium. 

Karrasch et al. provide evidence by DMRG data for 
the NN model that the equilibrium LL is universal also 
in nonequilibrium in the sense that Kq$ determines 
nonequilbrium power laws j^. The underlying idea is 
that the relevant energy cutoff of the RG flow is oc 1 /t so 
that for t — > oo the low-energy LL governs the physics. 

We prefer the view point that there arc two energy cut- 
offs 1/t and the excited energy AE/L because AE/L acts 
similar to a temperature smearing out smaller energy dif- 
ferences. Thus eventually max(l/i, AE/L) should stop 
the flow. In this case, the initial high energy LL param- 
eters relevant close to the FS flow very little if AE/L 
is large enough and the dynamics continues to be gov- 
erned by the FS exponents. For instance, the flow of 
K at Gs to its equilibrium value of 1 is logarithmically slow 
so that any substantial cutoff keeps it close to K ayF s- 
Since at present this issue cannot be decided further an- 
alytical and numerical work is called for: The regime 
1/t < AE/L for various values of AE/L would help to 
distinguish the two hypotheses. 

Our second key result stems from the comparison of 
the ID data with exact DMFT data for a Bethe lattice 
Q. The common lore would expect a crucially distinct 
behavior due to the differing dimensionality and the in- 
tegr ability of the ID Hubbard model [2^. The scattering 
in ID is strongly restricted due to momentum conserva- 
tion I17H191 l2ll. l22i while this conservation is irrelevant 



Indeed, the GS and the FS power laws differ significantly 



at internal vertices in DMFT [3l|. Yet Fig. [7] shows qual- 
itatively very similar results for larger values of U. While 
for U ;$ W differences prevail, the results for larger U, 
displaying coherent oscillations in An(t) agree up to 
the point that the minima appear at almost the same 
values. We conclude that these coherent oscillations are 
driven by local physics where the U term dominates. The 
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FIG. 5: (color online) Jump An(t) for the half-filled Hubbard 
model in ID (black, solid lines) and on the Bethe lattice (gray, 
dashed lines from Ref. ||). Dotted line: Second order result 
An = 1 - U 2 t 2 /A. 



right panel in Fig. [3] supports this view: The period of 
the oscillations approaches T = 2-k/U of local Rabi oscil- 
lation between a singly and a doubly occupied site. The 
lattice acts as a damping bath. To stress that this is not 
simply an effect of the leading order in t the second order 
result An(t) = 1 — U 2 n(l — n)t 2 (n filling factor of one 
spin species) is included in Fig. [7] 

Summarizing, we first showed that at short and moder- 
ate times the dynamics after interaction quenches in one 
dimension can be approximately described by a Luttinger 
liquid with parameters relevant around the initial states. 
For longer times, numerics is inconclusive, but we ex- 
pect that for substantial excitation energy the Luttinger 
parameters hardly flow so that the relevant exponents re- 
main close to the ones at short and moderate times: The 
equilibrium values are not reached. 

Our second key result are unexpected similarities in 
the noncquilibrium dynamics after interaction quenches 
in Hubbard models in one dimension and on the infinitely 
branched Bethe lattice. In spite of fundamental differ- 
ences in intcgrability and restriction in scattering both 
cases reveal damped coherent oscillations for larger val- 
ues of U with almost identical frequency. This points 
towards a physics dominated by local processes. 

Both these findings fundamentally enhance the un- 
derstanding of nonequilibrium physics in one dimension 
which is an intensely studied focus these days. Further 
analysis, e.g., at longer times, is called for, but very dif- 
ficult to reach to date. 
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SUPPLEMENTAL MATERIAL 

TECHNICALITIES 

We are interested in momentum distributions after in- 
teraction quenches. They can be computed by Fourier 
transformation from the one-particle equal time propa- 
gators 



G(f,t) = (FS|c(f,t) c t (0,t)|FS) 



(7) 



where the expectation value is taken with respect to the 
nonintcracting Fermi sea |FS) which represents the initial 
state before the quench. The r stands for a site on the 
lattice under study, here a chain in one dimension (ID). 
The time-dependent operators c and c 1 are represented 
by the following ansatz 



(8) 



where denote various superpositions of particle 

(hole) creation operators. For instance, the single parti- 
cle creation T> is given as 



(9) 



with time-dependent prcfactors ho(5,t). The possible 
shift S which occurs is given by the maximal velocity 
ti max in the sense of the Lieb- Robinson theorem [J, Q . 

To calculate the time dependence of the prefactors we 
use the Heisenberg equation dtA(r, t) = i H, A(f, t) for 

the time derivative of an operator A. On calculating the 
commutator H,c'(f,t) we encounter two cases. The 
commutation of the noninteracting part of the Hamil- 
tonian Hq leads to a shift of the fcrmionic operators 
in space, whereas the commutation with the interaction 
term H lnt may additionally create or annihilate particle- 
hole pairs Iterating this process then leads to 

the ansatz l[8]). With each commutation more terms with 
higher number of particles involved are created. Thus 
the amount of terms grows exponentially: At 11 loops 
we deal with up to 5 • 10 5 monomials in the Hubbard 
model and a set of differential equations with about 2- 10 7 
terms on the right hand side. In the spinless model with 
nearest-neighbor repulsion, we cope with up to 5 • 10 5 
monomials and a set of differential equations with 7 • 10 6 
terms on the right hand side. 

The differential equations of the prefactors can be 
solved numerically with the initial conditions ho(0, 0) = 1 
and hi(f, t) = Vi =/= 0. Because each commutation com- 
prises one order in time t a calculation with n commu- 
tations provides results for c)(t) which are exact up to 
order t n . Wc stress, however, that we do not use a scries 



expansion in t but solve the (truncated) set of differential 
equations numerically 

Due to the proliferating number of additional terms 
arising within the calculation we have to restrict ourselves 
to a finite number of commutations. As terms appearing 
during the last commutation for the first time lead to 
an overestimation of the weight loss for the one particle 
terms, we omit them to improve the convergence. A cal- 
culation with n commutations performed in this way is 
called an n-loop calculation. 



DYNAMICS IN THE HUBBARD CHAIN 
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FIG. 6: (color online) Solid lines: An(t) for the Hubbard 
model at quarter-filling for various values of U. Dashed lines: 
An(t) given by the ground state (GS) exponents. Dashed- 
dotted lines: An(t) calculated with the Fermi sea (FS) expo- 
nents. 



Figure [B] shows the jump An(t) of the momentum dis- 
tribution at the Fermi vector for the Hubbard model at 
quarter-filling. We compare the data for the microscopic 
model to the power laws for Luttinger liquids with two 
different sets of exponents: The ground state (GS) ex- 
ponents are taken from exact results; the Fermi sea (FS) 
exponents are taken from bosonization around the Fermi 
sea. As can be seen in the upper panel, the curves for the 
FS and the GS power law behave very similarly for small 
U ;$ 0.3IF. For larger values of U the time evolution of 
the jump is well described by the FS power law while the 
GS exponents do not fit the slope. 

At U = LOW the agreement worsens, which is ex- 
plained by the breakdown of the Luttinger liquid descrip- 
tion in terms of bosonic modes without interaction due 
to neglected scattering processes such as back scattering 
and umklapp scattering and due to the neglected curva- 
ture of the single-particle dispersion. 
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COMPARISON TO An(f) ON THE BETHE 
LATTICE 

Comparing the ID results with the results for the 
Bethe lattice with infinite coordination number obtained 
from dynamic mean field theory (DMFT) 0], the curves 
show a surprisingly good agreement for larger values of 
the interaction U. As can be seen from the upper panel 
of Fig. [71 the curves for small U only agree for very small 
times in the leading quadratic order in t. For U > W the 
DMFT and the ID curves agree qualitatively. Their min- 
ima appear even quantitatively at about the same times. 
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FIG. 7: (color online) Jump An(t) for the half-filled ID Hub- 
bard model for various values of the interaction U. Solid lines: 
An(t) for the ID Hubbard model in 11 loops. Dashed lines: 
DMFT data for the Bethe lattice taken from Ref. [j| 



